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Abstract. The Navier-Stokes-Voigt model of viscoelastic incompressible 
fluid has been recently proposed ffularization of the three-dimensional 

Navier-Stokes equations for the purpose of direct numerical simulations. Be- 
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sides the kinematic viscosity parameter, u > 0, this model possesses a reg- 

2 

ularizing parameter, a > 0, a given length scale parameter, so that ^ is 
the relaxation time of the viscoelastic fluid. In this work, we derive several 
statistical properties of the invariant measures associated with the solutions 
of the three-dimensional Navier-Stokes-Voigt equations. Moreover, we prove 
that, for fixed viscosity , z/ > 0, as the regularizing parameter a tends to 
zero, there exists a subsequence of probability invariant measures converg- 
ing, in a suitable sense, to a strong stationary statistical solution of the 
three-dimensional Navier-Stokes equations, which is a regularized version of 
the notion of stationary statistical solutions - a generalization of the con- 
cept of invariant measure introduced and investigated by Foias. This fact 
supports earlier numerical observations, and provides an additional evidence 
that, for small values of the regularization parameter a, the Navier-Stokes- 
Voigt model can indeed be considered as a model to study the statistical 
properties of the three-dimensional Navier-Stokes equations and turbulent 
flows via direct numerical simulations. 

MSG: 76D06,76D05,76F20,76F55,76A10 

Keywords: Navier-Stokes-Voigt, Navier- Stokes- Voight, inviscid regulariza- 
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1 Introduction 

In this work, we consider the three-dimensional Navier-Stokes-Voigt equa- 
tions (NSV), a viscoelastic incompressible fluid model introduced by Os- 
kolkov in [30], and which was proposed by the authors of [2] as a smooth 
inviscid regularization of the 3D Navier-Stokes equations (NSE) for the pur- 
pose of direct numerical simulations (DNS). More speciflcally, we consider the 
NSV model subject to periodic or no-slip boundary conditions, and driven 
by a given force fleld f. The velocity vector fleld, v(x, t), and the scalar 
kinematic pressure, p(x, t), are governed by the system of equations 

(9t(v - a^Av) - z/Av + V • Vv + Vp = f, x e n, 
V-v = 0, xGil, 
v(x,0) =vo(x), xGfi, 
v(x, t) =0 X G dQ, or v(x, t) is periodic; 
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in the smooth domain Q C M^, in the case of no-shp Dirichlet boundary 
condition, or with basic periodic domain Q = [0, L]^ C M^, when equipped 
with periodic boundary conditions. Here, a > is a given length scale 
parameter, and z/ > is a given kinematic viscosity, such that a^/z/ is the 
relaxation time of the viscoelastic fluid. 

In the periodic case, we assume that the driving force, and the initial 
velocity uq satisfy 



The Navier-Stokes-Voigt model of viscoelastic incompressible fluid, ([T]), 
(sometimes written as Navier-Stokes-Voight) was introduced by Oskolkov in 
[30j, and pointed out by O. Ladyzhenskaya as one of the reasonable mod- 
ifications of the Navier-Stokes equations, see [26]. In [30], A. P. Oskolkov 
studied and proved its solvability in different functional spaces. We remark 
that these equations behave like a damped hyperbolic (pseudo-parabolic) 
system, see [22], and, therefore, the solutions do not experience fast (in- 
stantaneous) smoothening of the initial data, as it is for parabolic systems 
like the Navier-Stokes equations. This fact could prevent the NSV model 
from being a reasonable modification of the Navier-Stokes equations, how- 
ever, since we are proposing it as a model for direct numerical simulations 
of turbulent flows in statistical equilibrium, i.e., after the solutions reach the 
global attractor, we are mainly interested in its long time behavior. Indeed, 
it was proved in [21], that solutions in the global attractor are smooth, if the 
forcing field is smooth enough, even for initial data satisfying only finite ki- 
netic energy and finite enstrophy (i.e. bounded in the Sobolev H^-norm). In 
particular, in [21], it is shown in the periodic case, that if the forcing field is 
analytic, then the global attractor consists of analytic functions. This result, 
in conjunction with results proved in this work, proves that if the forcing 
field is smooth enough, then averaged structure functions, with respect to an 
invariant measure for the NSV flow, display exponential decaying tail. 

As it was observed above, the NSV model presents an extra length scale 
associated with the viscoelastic properties of the fluid, the parameter a, 
besides the well known Kolmogorov length scale, rj, (see, e.g., [18] for its 




and it is easy to see that it implies 
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definition), which is usually associated with the smallest scales of motion in 
turbulent flows. In [22], based on numerical simulations of the Sabra Shell 
model, it is observed that for a large range of values of a > 0, i.e. a ^ rj, 
there are two distinct regions associated with the inertial range of the energy 
spectrum for the NSV model. The flrst one obeying the celebrated Kol- 
mogorov k~'^^^ power law (with anomalous correction), followed by a second 
range of wavenumbers, where energy condensates, and it is simply equipar- 
titioned. The range of the second power-law, however, slowly disappears as 
a/rj decreases, restoring the usual Navier-Stokes inertial range regime, when 

The numerical simulations of the Shell model presented in [23] also sug- 
gest that for a large range of values of a > 0, when a ^ t], small scales 
velocity fluctuations are vigorously damped, due to a slowness of the energy 
transfer timescales associated with these small length scales. This has an 
extraordinary effect in reducing the energy dissipation rate intermittency, 
characterized by violent fluctuations away from the average in the energy 
dissipation rate. It is therefore suggested in [23] that by tuning the param- 
eter a, one may attenuate the strong velocity fluctuations related to the 
intermittent events, reducing thus the stiffness of DNS of turbulent flows, 
with only a small effect on the energy containing scales. This reduction of 
stiffness in the NSV model has also been observed in numerical experiments 
in implementing the NSV model for inpainting [S|. 

The numerical experiments in [2S] were performed using the Sabra shell 
phenomenological model of turbulent flows. However, because these models 
do not present any spatial structure, it is not clear if such properties persist 
for full direct numerical simulations of 3D turbulent flows. In fact, a large 
class of viscoelastic flows present anomalous behavior even for very small 
viscoelastic parameters, see, e.g., [2S|, and such phenomenon could be present 
in flows governed by these equations. 

In this work, we rigorously establish several statistical properties of the 
energy spectrum of the 3D NSV model that were observed in the shell model 
simulations in [25], using the notion of invariant measures. For example, 
by considering the results of [21], we can rigorously justify the exponen- 
tially decaying dissipation range observed in the simulations presented in 
[25] . despite the fact that these equations behave like a damped hyperbolic 
(pseudo-parabolic) system [22] . 

Concerning our main result, we start by defining a notion of strong sta- 
tionary statistical solution of the Navier-Stokes equations, which is a reg- 
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ularized version of the concept of stationary statistical solution that was 
introduced and investigated by Foias in [TU], [H] (see also [IS])- This notion 
is a generalization of the concept of an invariant measure for the semigroup 
generated by the three-dimensional Navier-Stokes equations. 

Our main result states the following: Fix u > 0, and let /x"" be a given 
sequence of invariant measures for the semigroups generated by the NSV 
model. Then, there exists a subsequence, , weakly converging to a strong 
stationary statistical solution of the Navier-Stokes equations, as — >■ 0. 

This fact supports the observations reported in [25j , providing a rigorous 
evidence that, for small values of the regularization parameter a, the NSV 
model approximates several statistical properties of the 3D Navier-Stokes 
equations, and, therefore, may be used as a rehable subgrid scale model for 
direct numerical simulations of turbulent flows. 

2 Mathematical framework 

We will follow the standard functional formulation for the Navier-Stokes 
equations, see, e.g., [3|, [151 ES], and the functional formulation for the NSV 
equations used in [22]. We denote by V, in the no-slip case, the set 

v = {vec^{ny, v-<^ = o}. (2) 

In the periodic case, we define V as 

V = {All the three-dimensional vector valued trigonometric polynomials cp 
with basic periodic domain f2 = [0, L]^; V ■ v? = and J^ip dx = . 

(3) 

The two fundamental functional spaces in this work are defined by 

H = closure of V in {L'^{n)Y , 

and 

y = closure of V in {H\n)Y . (4) 
The inner products in H and V are denoted, respectively, by 

(u, v) = / u(x) ■ v(x) dx, ((u, v)) = / 
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and the associated norms by |u| = (u, u)^/^, ||u|| = ((u, u))^/^ (the latter is a 
norm thanks to the Poincare inequahty ([H]), below). 

We denote by Plh the (Leray-Helmholtz) orthogonal projector in L^(f2)^ 
onto the subspace H. The Stokes operator is defined by Am = — PlhAu, 
from D{A) = H"^ nV into H. The norm in D{A) is 

||u||^(^) = |Au|, VugD(A). 

Moreover, because the inverse of the Stokes operator is a positive self-adjoint 
compact operator in H, we can define its powers A'^, s G M, with domain 
D{A'). We have V = D{A^/^) and its dual V = D{A-^/^), see, e.g., [3j, [15], 
or [3S]. 

The term B{u, v) = Plh((u ■ is a bilinear term associated with the 

inertial term, which satisfies 

B{u,v)eV', u,^reV, (5) 

where V denotes the dual space of V, (see, e.g., [3], '35]). Thus, taking the 
duality action in ^ x V, between the bilinear term, i?(u, v) and a vector 
field w eV yields a trilinear term (see, e.g., [51 155]). 

6(U,V,W) = (P(u, v), w)y/y, 

which is defined for u, v,w G V, where (■, ■)yy denotes the duality action 
between the spaces V, and its dual space V. An important relation for the 
trilinear term is the orthogonality property, see, e.g., [51 155]. 

6(u,v,v) = 0, (6) 

for u, \ E V. It follows from this relation the anti-symmetry property (see, 

e.g., [SI [23), 

6(u, V, w) = -6(u, w, v), (7) 

for all u, V, w G V. 

Moreover, because the inverse of the Stokes operator is a positive self- 
adjoint compact operator in H, then there exists a complete orthonormal 
basis of H formed by eigenvectors, {wj}j, with associated eigenvalues, Awj = 
AjWj, satisfying < Ai < A2 < . . . < Aj ^ cxd, as j ^ cxd, see, e.g., [3[IT5|[35] 
for details. In this geometry and setting, the Poincare inequality holds 

|vr<ArM|vf , (8) 
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where Ai is the first eigenvalue of the Stokes operator. 

We will need the following inequalities, which are consequences of Holder 
inequality, Poincare inequality, and well-known inequalities of Ladyzhen- 
skaya, see, e.g., [3], [15], or [35] : 

|6(u,v,w)| < C|u| ||v|| |w|^/Vw|^/^ < C|uPvPw|, (9) 

for every u E H, and v,w G D{A). Also, 

|6(u,v,w)| < C|u|i/^|Au|^/^||v|| |w| < C|Aupv||w|, (10) 

for every u, v G D{A), and w E H. 

Throughout this work, if not otherwise stated, we assume a forcing field 
i E H. If we apply the Leray-Helmholtz operator, Plh, to the first equation 
of the system ([1]), we obtain the following equivalent functional differential 
equation 

^ + = F(v) := (f - uA^^ - i?(v, v)), v(0) = Vq. (11) 

Global existence and uniqueness of ffTTl) was first studied in [30] , where it was 
established that the system ([1]) generates a continuous semigroup, S"{t) : 
V ^ V, t E M"*", i.e., a one-parameter family of maps {5'"(t)}^>Q, satisfying 
the properties of a continuous semigroup, such that the solution of ([T]) (or 
equivalently, of (fTTl)). satisfies v(t) = S'(t)vo, for every vq G V. In |2J, global 
regularity of (ITT]) was also proved for the inviscid model, i.e. when u = 0. 
Moreover, higher-order regularity of the inviscid case, i.e. the Euler-Voigt 
model, is established in [2l]; furthermore, a new blowup criterion for the 
three-dimensional Euler equations, by means of this inviscid regularization 
is introduced. 



3 Energy budget for the NSV model 

The global regularity result for solutions of the Navier- Stokes- Voigt equations 
established in [30] implies that the following energy equality holds for every 
t G [0,oo): 

j^il H;t)f + y ||v(-,t)f ) = (f, v(-,t)) - u ||v(.,t)f , (12) 
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(in fact the energy equality and global well-posedness holds for all t G M). 
Similarly, the global regularity results established in [2] imply that the solu- 
tions of the NSV equations in the inviscid (i.e. Euler-Voigt model, u = 0) 
and unforced setting, f = 0, satisfy for every t G M: 

|(i|v°(..«)r + !^||v»(..«)||Vo. (13) 

Therefore, the conserved quantity in the inviscid and unforced setting of the 
NSV (i.e. Euler-Voigt) model is 

C ■'■ I |2 I II l|2 /-I.N 

So, = -M + — ||v|| , (14) 
which we call the a-energy. The quantity 

= ^ |V|' (15) 

is the usual kinetic energy. We remark that while the kinetic energy is for- 
mally conserved for the inviscid and unforced incompressible Navier- Stokes 
equations, this is not the case for the inviscid and unforced setting of the 
NSV model. Those differences in the conserved quantity for the inviscid 
and unforced case are the main source of deviation between the statististical 
properties of the NSV model and of the Navier-Stokes equations observed in 
|25j . In Section 6, we establish a result concerning the convergence of the 
averaged kinetic energy of the NSV model as the parameter a tends to zero. 

Let Wj be the orthonormal basis of H composed by the eigenf unctions of 
the Stokes operator, A. For a vector field v G -ff, we define the component 
Vk, for a wavenumber k, by 

We also define the component v^' k" by 

Vk',«;" = X Vk. (17) 

k'<k<k" 

Then, we can write the projected Navier-Stokes- Voigt equations on the shell 
[k', k") 

^(vk',k" + a^AvK',K") + lyAv^f^^'i + B{v, v)«;/_k" = i^^^^n. (18) 
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Let us now obtain the a-energy budget for the shell [k', n"). Multiplying 
in LP' by v^'^k"? "we obtain 

^^(|Vk',k"|^ + oi^ \\^k',k"\?) + ^ \\^k',k"\\^ = &(V, V, V«^,«//) ^^^^ 

- 6(v,v,v^,,^/) + (f,v«.,«") = [e",(v) - e"„(v)] + (f,v^/,«")> 

where 

e:(v) :=e:'-(v)-e:'-(v) 
is the net rate of a-energy transfer at k, and 

e"'^(v) := -(5(v«,,«,v«,,^),v«,oo) 

represents the net rate of a-energy from the lower modes to the higher modes, 
and 

e"'^(v) := -(fi(v^,oo,v«,oo),v«,,^) 

represents the net rate of ct-energy from the higher modes to the lower modes. 
Remark. Since we are dealing simultaneously with a > and a = 0, then 
for k" = oo, all the calculations above are formal. However, they can all be 
rigorously recovered for the NSV case, i.e. when a > 0, via the Galerkin 
approximation procedure. For the Navier-Stokes case, i.e. when a = 0, if 
k" = oo, we can only state equations f|T8l) and f|T9l) in terms of inequalities, 
due to a possible lack of regularity for the three-dimensional Navier-Stokes 
equations, see, e.g., [15]. 



4 Invariant measures and stationary statisti- 
cal solutions for the NSV model 

A rigorous mathematical framework for investigating the statistical proper- 
ties of turbulent fluid flows in statistical equilibrium was first considered by 
Hopf in In [IHIIII], Foias established the notion of stationary statisti- 
cal solutions of the Navier-Stokes equations, which is a generalization of the 
concept of invariant probability measures for the semigroup associated with 
the solutions of the equations of motion. This notion is important because 
the semigroup generated by the Navier-Stokes equations is not known to be 
well-posed in V, and, therefore, the notion of invariant measure in V is not 
well-posed either. A detailed discussion of this issue can be found in ^3]. 
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Because of the global regularity of solutions of the NSV model in the 
space V obtained in [30], one can refrain from using the abstract notion of 
stationary statistical solutions, and work only with the more familiar notion 
of invariant measure. 

However, because the main goal of this work is to show convergence results 
of the statistical properties of the NSV model to the corresponding ones of 
the Navier-Stokes equations, as a — > 0, we will show first that invariant 
measures of the NSV model are stationary statistical solutions of the NSV 
model, and using this fact, we prove a weak convergence theorem in the 
sense of stationary statistical solutions of the Navier-Stokes equations. This 
approximation result is the main reason why we are interested in deriving 
several statistical properties of the invariant measures associated with the 
NSV model in this section. 

We recall that a probability measure /x" on V is called an invariant mea- 
sure for the semigroup {>S'"(t)}^>g if 

Mi?) = /i((5"(t))-'i?), vt>o, 

for every Borel measurable set E in V. This definition is equivalent to say 
that 

[ $(u)(i/x"(u) = [ $(5"(t)u)d^"(u), V$ G L\fi). 
Jv Jv 

It is easy to see that the Dirac measure concentrated at the steady state 

solutions of the NSV model, which coincide with steady solutions of the 

Navier-Stokes equations, are invariant measures for the semigroup {5'"(t)}^>Q 

generated by the NSV equations. Therefore, because the set of steady states 

is nonempty (see, e.g., [3], [16] and [35]), we have that the set of invariant 

measures for {S°'(t)}^yQ is nonempty. 

We recall that the support of an invariant measure, /i", consists only 
of points that are nonwandering, i.e., for any u G supp/i", and any Borel 
set E containing u, there exists a sequence of positive times t„ — > cxd, such 
that S°'(tn)E n E ^ (/}, for all n G Z"*", see, e.g., [15j. This implies that the 
support of an invariant measure, /i", is included in the global attractor of 
the semigroup S°'(t). 

In [22], it was proved that S°'{t) : V ^ V has an absorbing ball in V, 
and it is an asymptotically compact semigroup, implying the existence of a 
global attractor in V. Moreover, the global compact attractor was shown to 
be bounded in D{A) and to have finite Hausdorff and fractal dimensions. 
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Furthermore, despite the fact that the NSV equations behave hke a 
damped hyperbohc (pseudo-parabohc) system, rather than a parabohc equa- 
tion, it was proved in [2T] that the 3D periodic NSV equations possess an 
asymptotic smoothing property. More specifically, let us define the Gevrey 
class of functions 

[ iez3\{o} 

(20) 

where Uj are the correspondent Fourier coefficients of the ^-periodic vector 
field u, r > 0, and r > 0. The space is equipped with the corresponding 
inner product 

(u,v),,. = (A^/V^''',A^/V^''')= urv,|jfV^I^'l<oo, (21) 

jez3\{o} 

for u, V G G^, and with the corresponding norm 

for u G One can prove that the space of real analytic functions C'^{VL) 
has the following characterization: 

T>0 

for any r > 0, see, e.g., [5], [27j for details. In [2]j|, it was proved that if the 
driving force, f, belongs to a Gevrey class of functions, G^^, for some tq > 0, 
then for every solution v(x, t) of the 3D periodic NSV equations, with initial 
data vq G V, there exists to > 0, and a function w(t) G L°°(0, oo, G^), for 
some r > depending only on |/|i,to) ^, and a, so that 

lim ||u(t) - w(t)|| = 0. 

t— >oo 

This result proves, in particular, that if f is analytic, then the global attractor 
consists of analytic functions. We recall that this technique of Gevrey class 
regularity was first introduced by Foias and Temam ^17] for proving the 
analyticity, in space and time, for strong solutions of the three-dimensional 
Navier-Stokes equations, for short time (see also for generalization of the 
technique to parabolic analytic equations). 

The following theorem is a straightforward application of the facts dis- 
cussed above: 
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Theorem 1 Let fi°' be an invariant measure for the semigroup {S°'{t)}^^Q 
generated by the 3D periodic NSV equations, ([T]) . // the forcing field f belongs 
to H , then the supp /x" is a subset of the global attractor, A, of the semigroup 
generated by the 3D NSV model, ([1]), which is bounded in D{A). Moreover, 
if i & Gl^, for some tq > 0, then supp fi" <Z A C. G^, for some r > 0, which 
implies 

Remark. We remark that the trigonometric polynomials, usually considered 
as forcing fields in the numerical investigations of turbulent flows, belong to 
the class of analytic functions. Therefore, our result applies to a wide set of 
simulations. In particular, it justifles the observed exponential tail in Sabra 
shell model simulations of the NSV equations in [25J (see similar results 
concerning the three-dimensional Navier-Stokes equations [H]). 

Now, we deflne the notion of stationary statistical solutions of the Navier- 
Stokes- Voigt equations inspired by Foias [IOl[TT], and [15]. 

Definition 1 A stationary statistical solution of the 3D Navier-Stokes- Voigt 
model, is a Borel probability measure /i" on V such that 

(2) [ ((/ + a^A)-^^'{u), f - uAu - B{u, u))c//i"(u) = 0, 
Jv 

for any test functional \I' G T", where is defined in Definition]^ 

(3) I [u \\uf - (f , u)) d^"(u) = 0, for every < Ei < 

J {Ei<\u\^ +a'2\\uf <E2} 
E2 < 00. 

Definition 2 We define the class of test functions to be the set of all 
real-valued functionals \E' = ^^(u) on H that are bounded on bounded subsets 
of V and such that the following conditions hold: 

1. For any u G D{A), the Frechet derivative \E''(u) taken in H along 
vectors in D[A) exists. More precisely, for each u G D{A), there exists 
an element in H denoted by \E''(u) such that 

' ^ ^ L^^^^O as |v| ^ 0,v G L'(A). (23) 
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2. \l/'(u) G H for all u G D{A), and the map u ^'(u) is continuous 
and bounded as a function from H into H . 

For example, we can take the cylindrical test functions : H —>■ W of the 
form ^(u) = ^/^ ((u, gi ),..., (u, gm)), where ip is a C^(IU"] R) scalar function 
on M"*, m G N, with compact support, and gi, . . . , gm belong to H. For this 
case we have 

ni 

^'(u) = djij{{u, gl), . . . , (U, gm))gj, 

i=i 

where djip denotes the derivative of ip with respect to the j-th variable. In 
this case, it follows that \Ef'(u) G H since it is a linear combination of the gj . 

The class of test functions are broader than the class considered in 
the usual definition of stationary statistical solutions of the Navier-Stokes 
equations, see, e.g., [Hj. This is made possible by the fact that solutions of 
the NSV model ([1]) are globally regular in V, see [2l 130] . 

Condition 1 in Definition [1] implies that the support of the stationary 
statistical solutions are included in D{A). Condition 3 is a local energy 
balance equation, which implies that the stationary statistical solutions have 
supports which are bounded in V, see Proposition [T] below. Moreover, in 
Corollary [1], we show that for fixed cto > 0, Condition 3 implies that the 
supports are uniformly bounded H, for all a G (0, ao]. 

We also remark that because we consider weakly converging subsequences 
of probability measures in the Section 6, we work with the weak topology 
of V. However, because is a separable Hilbert space, the Borel cr-algebra 
generated by the weakly open sets coincides with that for the open sets in the 
strong topology of V, see, e.g., |7]. Thus, we identify these two probability 
spaces in the rest of the work. 

Now, we prove a result concerning bounds for the supp /i° in the space 
V. This result will be important because we will derive from it a uniform 
bound, with respect to a, in the if-norm. 

Proposition 1 Let /i" be a stationary statistical solution of the Navier- 
Stokes- Voigt equations. Then, 

supp /i" C := I u G |u|^ + \\uf < (A^^ + a^)^- I . (24) 
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Proof. We follow the arguments used in [T3] for the NSE case. It follows 
from Definition [T] item (3), that if 

r := {u G V] El < \uf + ||uf < E^} , (25) 

then, by the Cauchy-Schwarz and Poincare, ([8]), inequalities, 



j ||uf ci/i"(u) < j |f| |u|d;u"(u) 

j ^ c//i°(u) j (^j ||uf rf/i"(u^ 



where Ai is the first eigenvalue of the operator A. Thanks to Poncare in- 
equality, (IH]), and Condition 1 in Definition [H the term in the right-hand side 
of the last inequahty is bounded, and therefore, 

|f|2 

- ^ 1 ^/^"H < 0- (26) 

Notice that if we choose, in the set F, Ei =0, and let E2 oo, and using 
the fact that /z" is a probability measure, we find that 

2 

/ ||uf rf/x"(u)<4L. (27) 
Jv ^ M 

Let < e < 1, we define 
and 

By Poincare inequality, ([8]), we have that |u|^ + ||u||^ < (A|f^ + «^) ||u"^ 
and, therefore. 



I |2 I 2 II l|2 |p|2 
, ,,2 U + a U I . 
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for all u G := {u G V\ E{ < |up + ||uf < E^]. We obtain, by apply- 
ing (EE]) to the set ^^ that /i°(r') = 0. Letting e i 0, we conclude (El). □ 

Now, we define a set that will be often used throughout the work 

:= |u G V- |u|' < (Ar' + . (28) 

The following corollary results immediately from Proposition [TJ 

Corollary 1 Let /i" he a stationary statistical solution of the Navier-Stokes- 
Voigt equations, then 

supp/i"c5„. (29) 

Next, we show that for every a > 0, invariant measures for the semigroup 
S'^it) : V ^ V generated by the NSV model, ([T]) (or flTTl) ). are stationary 
statistical solutions for the NSV model. The proof is in the same lines of a 
related result for the two-dimensional Navier-Stokes presented in |15j . 

Proposition 2 Let /i" he an invariant measure for the semigroup S°'{t) : 
V V generated hy the NSV model, (^!) (or ( ITTi) ). Then, /i° is a stationary 
statistical solution of the NSV model, in the sense of DefinitionUl 

Proof. Let fi" be an invariant measure for the semigroup S^lt) : V ^ V 
generated by the NSV model. Then, by Theorem [H supp /i° is included in 
the global attractor. A, of S"", which is bounded in D{A). Therefore, the 

1 1 2 

function u \\^\\d{A) bounded on the supp /i", which implies that 

/ l|u|lz)(A) '^/^"(^) ^ (^0) 

Therefore, Condition 1 in Definition [1] is satisfied. Now, we prove Condition 
3. Let < El < ^2 < oo, and let T = {ueV;Ei< \uf + ||uf < E2}. 
First, we observe that thanks to the Cauchy-Schwarz inequality, the Poincare 
inequality, ([8]), and to fl30|) . the map u — > (z/ ||u||^ — (f , u)) belongs to L^(/i"). 
Since /x" is an invariant measure, we have 

^(z/lluf -(f,u)) d^,-{u) = j^{u\\S"{t)nf ~ {f,S"{t)u)) rf/x-(u). 
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for all t > 0. Now, take the average with respect to t over [0,T], and use 
the fact that the left-hand side of the above equation is constant in time, to 
obtain 



)dt 



I (z/||uf -(f,u)) rf/x"(u) = l r I (H|5"(t)uf -(f,5°(t)u)) rf/x"(u 

= y"^(z/||^°(t)uf -(f,^"(t)u)) dtrf/x"(u), 

(31) 

where in the last step we are allowed to use the Fubini theorem, (see, e.g., 
[T5]). because the integrand is bounded in F, jointly continuous on [0, T] x V , 
and belongs to L^{\Q, T]xV), for all T > 0. By the energy equality f[T2]) . we 
obtain 



^(z/||uf -(f,u)) rf/i"(u) 

= ^ ^ (lu^ + ||uf - |S°(r)u|^ - a" ||S"(T)u||^) rf/x"(u) 



(32) 



for all T > 0. Since t S'"(t)u is bounded in V over [0, oo), then by letting 
T — > oo we obtain by the Lebesgue dominated convergence theorem that 



^(z^lluf -(f,u)) rf/x°(u)=0. 



which is Condition 3 in Definition [H For every m G Z"*", let us define the 
projector Pm : H ^ H hy 

m 

i=i 

where {w^} is an orthonormal basis of H consisting of eigenvectors of A. We 
now prove Condition 2 of Definition [1] Let $ G T° be a test functional. 
For m G Z^, we define $m(u) := ^(Pmu). It is easy to see that $m(u) = 
$'(Pmu)Pm is a functional on V. Because fi'^ is invariant, we have 

((/ + a2A)-iF(u),$:„(u))rf/i°(u) 

(33) 

((/ + a'A)-'F{S'^{t)u), $:„(5"(t)u))d/."(u) 
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Now, using the fact that the left-hand side is constant in time, we take the 
average with respect to t over [0,T], to obtain 



{{I + a'A)-'F{n),<!>'Ju))df,-{u) 

V 

T 



J^{{I + a'A)-'F{S-{t)u),<!>'^{S-{t)u))df,'^{n)dt (34) 
((/ + a^A)-^F{S''{t)u), $'^(5"(t)u)) d/i"(u). 



Jv 

T 



1 

'V ^ JO 

where we have used Fubini's theorem in the last step. 

Now, if u e D[A), then by the global regularity result of the NSV 
proved in [21 [30], see also [22], the function t — >• (J + a'^A)S°'{t)u belongs 
to C^((0, oo), H). Therefore, because the operator (/ + a^A)~^ : H ^ H is 
bounded, see Lemma [21 we obtain, by (fTTl) . that S"{t){u) satisfies 

^^"(t)(u) = (/ + a2A)-iF(5"(t)(u)) inH. (35) 

Therefore, 

|$^(5"(t)(u)) = {^'^{S-mn)), l^-WH) ^3g^ 
= i^'miS'-mn)), (/ + a'A)-'F{S-{t){u))), 

for all t > 0, and any u G D{A), where F is as in (fTTj) . Then, take the 
average with respect to t over [0,T], and substituting the result in the last 
line of f[MD. we conclude that 



((/ + a2A)-^F(u),$:,(u))rf/."(u)= 1 i[<|.™(5"(T)u)-$^(u)]d/i"(u). 

(37) 

Because $m is bounded, we obtain, as T — > cx), 

/ ((J + aM)-iF(u), $;„(u))rf/i°(u) = 0. (38) 
Jv 

Since = *^m(-PmAi)-Pm) have that 

/ (P„(/ + a^A)-'F{u), $'(P„u))d^°(u) = 0. (39) 
Jv 
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Now, because P^u ^ u, in if, as m oo, for all u G we have 

(P^(/ + a^A)-^F{u), ^'{Pmii)) ((/ + aM)-iF(u), $'(u)), as m ^ oo. 

(40) 

Therefore, because $' is bounded, we have by ([9]) that there exists a constant 
c > 0, which is independent of a, such that 

|(P„(/ + a2A)-iF(u),$'(P„u))| < c(||u||^(^) + |f|), 

where the right-hand side belongs to L^^djjf^), by (!30|) . Then, by the Lebesgue 
dominated convergence theorem, we obtain 

f ((/ + aM)-iF(u),$'(u))rf;u"(u) = 

, (41) 
lim / (P„(/ + a2A)-iF(u), $'(P„u))c//i"(u) = 0. 

Therefore, Condition 2 holds, completing the proof of the theorem. □ 

Corollary 2 The set of stationary statistical solutions of the NSV equations 
is nonempty. 

Proof. As we have observed earlier, the set of steady state solutions of the 
NSV, or equivalently of the NSE, is non-empty, see, e.g., jSj 135j- There- 
fore, the Dirac delta measures concentrated on the steady state solutions 
and all their convex combinations are invariant measures. Consequently, by 
Proposition [21 they are stationary statistical solutions of the NSV model. 

Remark. One could also prove the converse of Proposition [2], i.e., that every 
stationary statistical solution of the NSV model is an invariant measure, by 
following the arguments of Foias in (TUl [H] for the 2D periodic case, see also 
|15j . Only a slight modification in the treatment of the nonlinear term is 
necessary, but which can be circumvented by using that supp /x" is included 
in D{A), and is bounded in V . We do not present it here, since the semigroup 
generated by the NSV model is well defined in V , and therefore, we can use 
the more natural notion of invariant measures in our approximation results 
in Section 6. 

In the sequel, we use the symbol (. . to denote average with respect to 
an invariant measure of the NSV model i.e. 

(<^)°:= / <^(u)rf/i"(u); Vv9GLi(rf/i-). 
Jv 
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We can use the stationary statistical solutions formalism to show that 
the stationary Reynolds averaged equations hold for the NSV equations in 
this framework. Indeed, we can define the mean velocity field (u)" G D{A), 
thanks to Condition 1 in Definition ([1]), as follows 

((u)",v) = /(u,v)rf/i"(u) = / (u,v)d/x"(u), VvGZ}(A)'. 
Jv Jd{A) 

Moreover, thanks to (ITOl) and Condition 1 of Definition [3l we may define the 
average (-B(u, u))" & H hj 

((fi(u,u))",v) = /"(5(u,u),v)c//i"(u) = [ (fi(u,u),v)rf/i°(u), Vv G 
Jv Jd{A) 

The mean flow (u)° is a vector field on Vt with (u)" G D{A), while 
(i?(u, u))° G -ff. The next proposition shows that, assuming statistical equi- 
librium, the vector field (u)" also satisfies the Reynolds averaged equations 
in H. This should be compared to the NSE case, where these equations are 
known to be valid only in V' , (see, e.g., [T5]). 

Proposition 3 Let /x" he an invariant measure of the NSV equations. Then, 
the following functional form of the Reynolds averaged equations holds in H: 

z/A(u)° + (S(u,u))" = f. (42) 

Proof. We will use here the notation of Section 3, and set k\ = \i. Let ■?/' 
be a real-valued function with compact support on R. For any v G if, 
any finite wavenumber k > ni, and every u G D{A), the real- valued function 

$(u) = ^{{{I + a'Ay/'u, (/ + a'Af\^,,,)) (43) 

is a cylindrical test function, and it satisfies 

($'(u), 0) = ^'(((/ + a'Ay/'u, (/ + a'Ay/\^,,MiJ + v.,,.). 

for every (f) G D{A). Thus, considering (/) = (/ + a'^A)^^F{u.), which belongs 
to D{A), we obtain 

($'(u), {lWA)-'F{n)) = ^p'{{{IWAY/'u, (/WA)1/V«„«))(F(u), v«,,.). 
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By Proposition [21 /i" is also a stationary statistical solutions of the NSV 
model, therefore, by substituting the expression above in Condition 2 of 
Definition [H, we obtain 

[ + a^Afl^MXl + aM) V2v,,,.)) {(f , V,,,.) 

Jv 

-Z/(y4u, Vki,k) - 6(u, U, V«;^,k)} d/i"(u) = 0. 

Now, by virtue of Proposition [1], supp /i" is bounded in V , and therefore 
the map u G supp ^ ((/ + a^A)-'^/^u, (/ + a^74)^/^VKj_K) is uniformly 
bounded in supp yu", let us say by M > 0. Therefore, we can choose above, 
such that ^' = 1 at the interval [-M, M], therefore, ^/''(((/ + aM)^/^u, (J + 
Q;^y4)^/^VK^ k)) = 1; for every u G supp /i". This yields 

{(f , v«j,«,) - z/(ylu, v^j,«,) - 6(u, u, v«,^,«,)} rf/i"(u) = 0. 

For each fixed v G -ff, we observe that since supp C D{A), we have by 
the Cauchy-Schwarz inequality, and by ([9]), that 

|(f, Vki,k) - z/(v4u, Vki,«;) - 6(u,u, Vki,k)| < |f||v| + z^|v4u||v| + C|v4u|Vl- 

Thanks to Condition 1 of Definition [H we can apply the Lebesgue domi- 
nated convergence Theorem, and let ^ oo to obtain 

/ {(f , v) - z/(Au, v) - 6(u, u, v)} rf^"(u) = 0, 
Jv 

which yields the result. □ 



V 



5 Averaged Energy Budget for the NSV model 

In this section, we follow [15] to investigate the energy distribution scale-by- 
scale for the 3D NSV equations. For the Navier-Stokes equations, assuming 
that there exists an inertial range, i.e. an extensive range of wavenumbers 
where the viscous dissipation does not play a significant role, one can show 
that the energy simply cascades through these wavenumbers with a rate equal 
to the mean energy dissipation rate, e = i/(||u||^)". For the NSV equations, 
a similar scenario holds for the a-energy 

C ^ I |2 , "^1 ||2 

Cq, = - u H u . 

2 ' ' 2 " " 
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As it is usual in the studies of homogeneous turbulence, and using the 
notation of Section 3, we will consider the forcing f with finite eigenmodes, 

f= J2 (44) 
for the orthonormal basis, {wj}, of H consisting of the eigenvectors of A, 

with K>Ki = 

Now, we argue as in the proof of Proposition [2] to obtain an averaged 
energy balance equation. We skip the details, but we sketch the derivation: 
Using the fact that /i" is an invariant measure, using equation (fT8l) . averaging 
with respect to t over [0,T], and using Fubini's Theorem, we have 

^ (z^ l|u«',«"f - (f,u«'A") - K>iu) - e^„(u)]) dfi°{u) 

(45) 

where F is as in the proof of Proposition [21 Now, let n" > n' > k, and letting 
T — > oo in the above equation, we obtain the following balance equation 

u{\\u^',AY = {e:>in)-e:.iu)r. (46) 

The expression on the right-hand side of the last equality is the mean 
net a-energy transfer in the energy shell [k',k"]. Because the supp /i" in 
included in -0(^4) and is bounded in V, we can use the Lebesgue dominated 
convergence Theorem, and let k" — oo to obtain 

^(||u..,oof)" = (e.'(u))". (47) 

This expression shows that the net a-energy transfer is positive for every 
K, > K. We remark that a similar expression holds for the three-dimensional 
NSE case substituting equality by inequality, see, e.g., [15]. Therefore, we 
have proved 

Proposition 4 Let yu" be an invariant measure for the semigroup, {^'"(t)}^^^, 
generated by the NSV model, and let R < k' < k," . Then 

= i{e.'iu)r - (e."(u))")- (48) 
Moreover, for all n> K, we have 

y{w,oofr = M^)r. (49) 
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The result above shows that the mean net a-energy transfer is positive, which 
is an important consistency check for turbulence models. Moreover, if we 
assume that there exists a range of wavenumbers, k"], where the viscous 
dissipation term satisfies z/(||uk',k"||^)" -C 1, it follows from (BHl) that the 
a-energy transfer is nearly constant within this range, that is, there exists 
an inertial range for the NSV model. 

6 Statistical approximations of the 3D Navier- 
Stokes equations 

In this section we prove that given a sequence of invariant measures for 
the NSV equations, {/i""}, with the regularizing parameter q;„ converging 
to 0, as n — > cxD, there exists a subsequence converging weakly to a strong 
stationary statistical solutions of the 3D Navier-Stokes equations, a notion 
that we introduce in Definition 3 below. This is a slight modification of the 
notion of stationary statistical solutions of the 3D NSE introduced by Foias 
in [iniin], see also [15]. 

The main difficulty in proving this result is that the supports of the 
invariant measures, /x", of the NSV are not uniformly bounded, with respect 
to the parameter a, in the norm of V, as a tends to zero. This means that we 
cannot use in a straightforward manner the standard Prokhorov's Theorem, 
see, e.g. [M| Thm. 1.12]. This is because the whole space V is not metrizable 
when endowed with its weak topology. Therefore, we have to consider more 
general theorems concerning convergence of measures in arbitrary Hausdorff 
spaces. Fortunately, the family of measures considered in this section satisfy 
the tightness condition of and the weak convergence holds for our case. 

We start by defining the terminology, from [SB], that is relevant to the 
results presented in this section. A paving is a non-empty set consisting of 
subsets of a given set X. Af <Z2-^ is said to be a (0, U/, flc)— paving of X if 
N is closed under finite unions and countable intersections, and if in addition 
e A/". Similarly, A/" C 2^ is said to be a (0, U/, n/) -paving of X if A/" is 
closed under finite unions and finite intersections, and G A/". A paving /C is 
called compact [semicompact] if every family [every countable family] of sets 
in /C, which has the finite intersection property (the intersection of any finite 
subset of /C is non-empty) has a non-empty intersection (the intersection of 
all elements of K, is non-empty). A paving C separates the sets in K, if to any 



22 



pair Ki, K2 of disjoint elements in /C, we can find a pair Gi, G2 of disjoint 
sets in C such that Ki C Gi and K2 C G2- 

The results stated in [36j consider a Hausdorff topological space X, and a 
net, {x/3}, in X. The notion of net can be found in [36j, but it is not important 
for our present work because the space \^ is a separable topological space with 
its weak topology, and thus all statements can be done in terms of sequences. 
However, we will state a theorem appearing in fSB] in its original formulation, 
in terms of nets. We only recall that a net {x 13) 13^0 on a topological space X 
is compact if every subnet has a further subnet that converges. 

Let X be a Hausdorff space, and let £ and /C be pavings in X. We denote 
hy B = B{X,C), the smallest cr-field containing every set E G X for which 
L n E E C, WL G C We denote by A4+{X, C) the set of finite, non- negative 
measures defined on B{X,C). £, /C) denotes the set of measures in 

A1+(X, £) which are regular with respect to /C, i.e., /i G A^+(X, £, /C) if and 
only if n satisfies 

^i{E) = sup il^iK); K e IC, K C E} , \/E(lB. (50) 

From now on, we will consider the following axioms 

1. X is a Hausdorff topological space. 

2. /C is a (0, U/, nc) -paving of X. 

3. £ is a (0, U/, n/) -paving of X. 

4. K\G e IC,yK e IC,yG e C. 

5. C separates the sets in /C. 

6. /C is semi compact. 

Let fi E 7V1+(X, £), and let {jJ^^} be a net in Ai+{X, C). Then, we say 
that {jJ'^} converges to /i in the weak topology if and only if — > 
and liminf > fi{G), VG G £. 

Now, we are going to state the abstract result that we need to continue 
with our investigation. It is just the first half of Corollary 1 in [55] . 

Theorem 2 (Corollary 1, 136^ ) Assume that axioms I — VI are satisfied and 
let {fi^)i3^D be a net on £,/C) such that limsup^g^ /U^(X) < oc. If 

the tightness condition 

inf lim sup /(X \K) = Q (51) 
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holds, then {fi'^)/3(^D is compact. 

We will now translate all the preceding abstract definitions to our specific 
context. We consider the topological space X = V, where V is defined in @, 
endowed with its weak topology. Thanks to the Hahn-Banach Theorem, the 
space \^ is a Hausdorff topological space with the weak topology. Moreover, 
V is also separable with the weak topology, see, e.g., |Kj. This implies that 
there is no need to use nets in the application of Theorem [2] in this work, 
and we will consider only sequences in the subsequent results. 

We denote by /C(V), J^{V), C{V) and B{V) the pavings of V of weakly 
compact sets, weakly closed sets, weakly open sets, and weakly Borel sets, 
respectively. By Ai+iV, C(y)) we denote the space of finite non-negative 
measures defined on B = B(y, C(y)). It is easy to prove that the axioms 
I -VI are vahd for the above choice of IC{V), J^{V), C{V) and B{V). We 
omit the proof of these facts. 

The weak topology on A4+{V,B{V)) is defined via the paving C{V), and 
it is exactly the familiar weak topology in the well-known Polish spaces, i.e. 
fi^ converges in the weak topology to /i if and only if 



for every h : V M bounded continuous function. We also recall that the 
subset of Borel probability measures Aii is closed in /iA^{V) under the weak 
topology, (see, e.g., [7j), which implies that weakly converging sequences in 
Ail converge to elements in TWi. 

Now we will prove that every Borel finite measure defined on V endowed 
with the weak topology satisfies the tightness condition of Theorem [21 

Lemma 1 Let fi e M+{V,C{V)), then /i G M+{V, C{V),IC{V)). 

Proof. First, we notice that because V with the strong topology is a sepa- 
rable Hilbert space, the Borel cx-algebra associated with the strong topology 
is the same as the Borel cx-algebra associated with the weak topology, (see. 



Moreover, because V with the strong topology is separable, any strongly 
Borel probability measure, fi, is regular in the sense that for every strongly 
Borel set E, 




e.g., m)- 



li{E) = sup {fi{K); K is strongly compact , K C E} 



(52) 
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see, e.g., [7]. Now, because strongly compact sets in V are also weakly 
compact, we have that if is a weakly Borel set in V (which also implies 
that it is a strongly Borel set), then /i satisfies 

> sup {fi{K); K is weakly compact , K C E} 

> sup {fi{K); K is strongly compact , K C E} = fi{E). 

Therefore, 

fi{E) = sup {fi{K); K is weakly compact , K C E} , (54) 

which shows that /i G M+{V, C{V), /C(V)). □ 

Now, it is easy to see that Theorem [2], Lemma [H and the preceding 
discussions imply the following theorem: 

Theorem 3 Let K, denote the paving of V formed by the weakly compact 
sets. Consider the space V endowed with its weak topology, and let {fin} be 
a sequence of Borel probability measures in V . If the tightness condition 

inf lim sup fi^iV \ K) = (55) 

holds, then there exists a subsequence, also denoted by {/Xn}, such that /i„ ^ 
weakly m M + {V,C{V)). 

Because of the convergence result of our main theorem in this section 
occurs in the space V , endowed with its weak topology, the limit measure 
thus obtained will be defined on this space as well. Therefore, we need to 
define the following notion 

Definition 3 A strong stationary statistical solution of the 3D Navier-Stokes 
equations is a Borel probability measure fi on V such that 

(1) I ||u||^ (i/i(u) < oo; 
Jv 

(2) I (^'(u),f-Z/Au-fi(u,u))yy/d^(u) = 0, 

Jv ' 

for any test functional \1' G T, where T is defined in Definition^ 

(3) I [u \\uf - (f , u)] d/i(u) < 0, for every < Ei < E2 < oo. 

J Ei<\uf<E2 
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As remarked before, we may use the weakly Borel a-algebra in Definition 
[3] because is a separable Hilbert space, and therefore, weakly Borel sets 
coincide with the strong Borel sets, see, e.g., 

Definition 4 The class of test functions T is the set of functions : V M. 
of the form 

vl>(u) :=^((u,gi),...,(u,g^)), (56) 

where the function ip is a scalar valued function defined on M™, with 
m G Z"*", and gi, . . . , gm belong to V . 

Remark 1: We remark again that for f G i/, there exist steady state so- 
lutions of the 3D Navier-Stokes equations which belongs to D{A), see, e.g., 
[31 [35]. Therefore, every Dirac measure concentrated on a steady solution of 
the 3D Navier-Stokes equations which belongs to D(A) is a strong station- 
ary statistical solution. Moreover, if u(t) is a strong solution of the NSE, 
bounded in V for every t > 0, then the trajectory of u(t) is included in a 
closed ball in V , which is compact in the weak topology of V . Therefore, 
one can also generate a strong stationary statistical solution by the Krylov- 
Bogolyubov procedure for generating invariant measures; see, e.g., [15] 
for details on how to use this procedure to generate time-average measures 
and the Banach limit. 

Remark 2: We remark that strong stationary statistical solutions are de- 
fined over the space V with its corresponding Borel cr-algebra, while that 
the usual stationary statistical solutions, defined by Foias in [10], [11], are 
defined over the space H with its Borel cr-algebra. This justifies the denom- 
ination strong stationary statistical solutions. Indeed, because Leray-Hopf 
weak solutions of the 3D NSE are uniformly bounded in the if-norm, with 
respect to time, one can generate, by the Krylov-Bogolyubov procedure, a 
Borel probability measure defined over H , which is a stationary statistical 
solution, (see, e.g., However, because strong stationary statistical so- 

lutions are defined over V , the Krylov-Bogolyubov procedure can generate a 
strong stationary statistical solution only if applied to strong solutions of the 
3-D NSE, that are uniformly bounded, with respect to time, in the V^-norm, 
as described in the former remark. 

We now proceed with the computation of for test functions \1' G T. It 
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is easy to see that for \1/ G T, as defined in (jSHD, we have 

m 

Vu^(u) ■ = ^ dj%lj ((u, gi), . . . , (u, gj) (0, gj). 

i=i 

It is easy to see that — > Vu^&(u)-0 is a bounded hnear continuous functional 
in H, hence, by the Riesz representation theorem, there exists an element 
e H such that 

Vu^(u)-0= (^'(u),0), V0Gi/. 

This is the identification that is used in Definition [31 More precisely, 

m 

*'(u) = dji' ((u, gi), . . . , (u, g„)) g,. (57) 
i=i 

We now establish a lemma that will play a key role in the rest of this 
work. 

Lemma 2 For every a > 0, the operators (J + a^A)^^ : H ^ H satisfy 

\\{I + a'Ar'\\^^^<l. (58) 
Moreover, there exists a constant C > 0, independent of a, such that 

I ((/ + a^A)-^ - /)0| < Ca^Xl-'^lAyi (59) 
for every G DlA"^), and rj > 7/4. 

Proof. Indeed, let {wj} be a complete orthonormal basis of H formed by 
the eigenvectors of the operator A, and let < Ai < A2 . . . < Aj — > 00, as 
j ^ 00, be the corresponding eigenvalues of the operator A associated with 
the eigenvectors {wj}. Then, by the Parseval identity, for every a > 0, 



1 + a^Xj 



l(0,w,)|'< 



which proves that ||(/ + a'^A) ^Wjj^jj < 1. Now, for every G if, we have 
by the Parseval identity that 



((/+aMr'-/M^f;(^-i) 

j=i ^ J / 



i=i 
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Therefore, because Wj 



a: 



-2r?| 



\{Ay,wj)\'' < Xj^^'lAyi'', for all j > 1, 



we once again have by the Parseval identity that 



\{{I + a^A)-' -I)4,\ 



E 



1 + a'^Xi 



E 



1 + a^Xj 



2 

,2(r,-l) , 



(60) 



2(r,-l) ' 



The asymptotic behavior of the eigenvalues of the Stokes operator defined in 
Vl C M" is known to satisfy the Weyl-type formula; see, e.g., [T], [25] : 



A, 



(2vr) 



u;„,(r7, — 



2/n 



/c^/", as ^ oo. 



(61) 



where is the n-dimensional Lebesgue measure of fi, and ct;„ is the volume 

of the unit ball in M", which implies that the infinite series X^^li ( ^ ) 
converges for rj > 7/4, if n = 3. Therefore, by (j6Til . it is easy to see that 

|((/ + a^A)-^ - I)(f)\ < CXf^-''^a^\Ay\, 

Where = Er^i (t) • ° 

The following lemma, which is a trivial consequence of the compact em- 
bedding V GG H, will play a key role in the sequel of this work. We prove 
it here for the sake of completeness. 



u| is weakly continuous with respect to the 
' , where Ba is as defined in 



Lemma 3 The function u - 
V -topology, and hounded over Ba C V , a E 

Proof. The bound follows immediatelly from the definition of B^ in 
Now, we prove that u — > |u| is weakly continuous with the l^-topology. 
Indeed, suppose that uj u weakly in V, as j — > oo, but Uj does not 
converge strongly to u in H. Then, there exists e > 0, and a subsequence 
{uj^} of {uj} satisfying 



u 



u| > e, (62) 

for all k G Z+. Now, because u^^ ^ u weakly in V, as k oo, we have, by 
the compact embedding V GG H, that there exists a further subsequence 
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|ujj.^| of {uj^} converging strongly (of course, also weakly) in to a vector 

field \ E V, a.s i ^ oo. Because weak convergence in V implies weak 
convergence in H, u^^^ ^ u weakly in if , as £ oo. By the uniqueness 
of the weak limit, u = v. Therefore, |uj^^ ~u| ^0, as£-^oo, which 
contradicts flB2]) . This proves the lemma. □ 

In order to prove the convergence results in the sequel of this work, we 
need first to consider a subclass of smoother test functions, 7i C T, consisting 
of functionals satisfying fl56l) . so that the vector fields gj appearing in their 
definitions possess a higher regularity, namely, gj G V, where V is given in 
© or da]). 

Lemma 4 Let ^ G Tl . Let ao > be fixed, and let B^o be as defined in (l28l) . 
Then \l''(u) G 0,^^^+ ^i^^)' ^'^^ there exists a constant, C > 0, depending 
only on £, t] and a^, such that 

\A^—m'{v.)\<C (63) 

holds for all ue B^^, £e Z+, j = 1, 2, 3, and r]>0. Consider F^"" : V ^ R, 
i = 1, 2, 3, a > defined by 

Ff(u) = ((/ + a^A)-V(u),f), (64) 



F,"{n) := (((/ + a'A)-'^'{n), u)) = {A{I + a'A)-'^'{n), u), (65) 

and 

F^iu) = (Bin, u), (/ + a'A)-'^'in)))v',v. (66) 

These three maps are well defined for u G B^, weakly continuous in V, and 
bounded in B^, with uniform bounds for a G [0, ao]. 

Proof. We first check the bound fl63|) . Since ip is of class C^, there exists 
C > such that for all u G B^a, we have 

|9,V'((u,gi),...,(u,gJ)| <a (67) 

Moreover, because for every vector field g G V, we have that 



g G fl D{A^), (68) 
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and l^^^gl < C, for some C > 0, for every ^ G Z"*", and j = 1,2,3. 

Therefore, we obtain fl63|) by simple inspection of fl57|) . fl671) and fl68|) . 

Now, we will prove that the maps i^"(u) are weakly continuous in V, and 
uniformly bounded with respect to a G [0, Oq] and u e -Bog. 

Throughout this proof, let {u^} C denote a sequence converging 
weakly in to a vector field u in V. Of course, this also implies that 
converges weakly in H to u. By now, it is easy to see by inspection of (l57|l . 
and by using (!58l) and ( |67l) . that the weak convergence in H of Uj to u implies 
that 

(/ + a''A)-^-^AW(uA ^ (/ + aM)-i^A''^'(u) (69) 

dxk dxk 

strongly in if, as j — > oo, for every k G Z+, and 77 > 0. Now, because the 
sequence of real numbers given by the inner product in if of a strongly con- 
vergent sequence in H, and a weakly convergent sequence in H is convergent, 
we have 

FfK) = {{I + a'A)-'^\u,),f) ^ ((i + a2A)-iM/'(u),f) =Ff(u), (70) 
and 

F,"in,) = ((i + a'A)''A^'{n,), u,) ^ ((i + a2A)-iAvl/'(u), u) = ^^-(u), 

(71) 

as j —>■ 00. This shows that and F2 are weakly continuous in V. The 
uniform bounds follow by ( 163|) . and by inspection of ( 164|) . 

Now, concerning the map F^, it is easy to see by the strong convergence 
in H of Uj to u established in Lemma [3l and by the strong convergence in 
H stated in ([69]), that 

5(u„ (i + a2A)-i^'(u,)) ^ fi(u, (i + a2A)-i^'(u)) (72) 

strongly in H, a.s j —>■ 00. Hence, by the anti-symmetry of the trilinear term, 
([7j), by the strong convergence of uj to u, and by (!72|) . it is easy to see that 



F3"(u,) = -(u,, B{u„ (i + a'Ar'^'iu,))) 
^-(u,i?(u,(i + a2A)-V(u)))=F3"(u) ^ ^ 



as j —>■ 00. This shows that F^ is weakly continuous in V. The uniform 
bounds follow from the following inequality, which is an easy consequence of 

0, 

|6(u,u,v)| < C|uHVv|^^, (74) 
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for every u ^ V, and v G and from the following Sobolev's imbed- 

ding inequality in 3D: 

\if\^^ <C\A<f\, ifeD{A). (75) 

Indeed, because dx^'^'iu) G D{A), for all u G -Bap, and j G {1, 2, 3} we have 
by ([63D, and by ([MD, that 

|F3"(u) I < C lul^l (/ + a'A)-'Vx^'{u)\^^ 

= C sup |(/ + a2A)"i^^'(u)Uoo|u|' 
je{i,2,3} <^^j 

<C sup |(/ + a'A)-M^^'(u)| |u|' "^^^^ 

je{l,2,3} OXj 

d ^, , ,2 / ^ I ,2 



<C sup |A— ^'(u)| |ur < C„Ju 

ie{i,2,3} 



where Cqq is uniform for a G [0,a;o]. This implies the uniform boundedness 
for for u G 5^,0 . □ 

Lemma 5 Lei ^ G T^. Consider Fi defined as in fl64l) . (!65|) and (!66|) wi/i 
this choice of'^. Then, for any sequence {/i""} of probability invariant mea- 
sures of the NSV, with a„ ^ 0, as n ^ oo, we have 



{Ft"{u)-F^iu))df,-"{u) 



0, as an 0, fori = 1,2,3. (77) 



Proof. Fix ao > 0, and now notice that because /i"" are stationary statis- 
tical solutions of the NSV model, by Proposition [2|, we have by ([1]) that the 
supports of are included in Ba^, for every a„ < ao. From now on, we 
restrict ourselves to < < ao- 

Now, by inspection of by (!59|) . and by (!63|) . it is easy to see that for 
every u G -Bq-q, there exists Cq,,, > 0, such that 

|F-(u) -Fi^u)! < C|f| |((/ + a^A)'i -/)vl/'(u)| 
<C|((/ + a^A)-i-/)M/'(u)| <CAr^a^|A''M/'(u)|<C„,a^. ' 

Now, by inspection of (|^ . by (jSHD, and by (1^ . it is easy to see that for 
every u G -Bq-q, there exists Cq,,, > 0, such that 

|F-(u) - FO(u)| < C |u| |((/ + alA)-' - /)Avl/'(u)| 
<C,J((/ + a^A)-i-/)AvI/V)| <CA}-^a^|^,+i^.(^)|<^^^^^ I J 
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For F^, because "^'{u) G D{A), for all u G Bag, and j G {1, 2, 3}, we have 
by dZD, dSlD, dMl) and dZg) that for a„ G [0,ao], 

|F3-(u) - F^{u)\ < |(u, B{u, ((/ + - lW{u)))\ 

<C\uf\{iI + alA)-'-I)V.^'{u)\^^ 
=<C\u\' sup |((/ + a2A)-i-/)a..M/'(u)L 



je{i,2,3} 



<c. 



ao sup 

ie{i,2,3} 



d 



{iI + alA)-'-I)A—^'{n) 



m 



< CaoalXl sup 

i6{l,2,3} 



Therefore, by fl78|) . fl79|) . flHOj) . and since /i'' are probability measures, we have 



/ |i^""(u)-i^°(u)M/x-"(u)<C,,a^ 



(81) 



The result follows as n oo. □ 

Now, we prove the main result of this work. 

Theorem 4 Given a sequence of invariant probability measures of the 3D 
Navier- Stokes- Voigt model, {/i""}, with an 0, as n oo, there exists a 
subsequence, denoted also by {/i""}, and a B or el probability measure fi on V, 
such that 

lim / $(u)rf/i""(u) = / $(u)rf/i(u), (82) 

Jv 

for all ^ : V ^ M. weakly continuous bounded real-valued functionals. Fur- 
thermore, the weak limit measure fi is a strong stationary statistical solution 
of the 3D Navier-Stokes equations. 

Proof. Because {/i""} are invariant measures, by Proposition [H we know 
that they are also stationary statistical solutions of the NSV model, and, 
therefore, they satisfy Conditions 1, 2 and 3 in Definition [H We prove first 
that the sequence {yu""} satisfies the tightness condition (1551) . Rewriting 
more explicitly inequality (l26l) . we obtain 



|uf d/i""(u) < 



Ei<\u\^+al\\u\\^<E2 



(83) 
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for every Ei,E2 > 0. Now, fixed i? > 0, choose e[^^ = Q^B?, and let 
E2 ^ oo, which is possible because supp/i"" is bounded in V. This yields 

[ ||uf rf/i""(u) < /" ||uf rf/i""(u) < (84) 

jR2<\\uf J {alR'^<\u\^+al\\uf} ^ M 

The first inequality above holds because of the relation 

{u G V; \\uf > R^} C {ueV; \uf + al ||uf > a^R^} , 

1 1 2 m 

and the fact that ||u|| is positive. Therefore, denoting by Bji2 the closed 
ball with radius R^ in V, we have by flM|) that 



1 |f|^ 
R^^i 



Because the closed balls Bfi2 are weakly compact in V, this yields the tight- 
ness condition flS^ . necessary for Theorem [31 Therefore, it is easy to see 
that Theorem [3] applies to the sequence , yielding the existence of a sub- 
sequence, denoted also by {/i""}, and a Borel probability measure fi on V, 
such that ([82]) holds. 

Now, we prove that is indeed a strong stationary statistical solution 
of the NSE. Condition 1 of Definition [3] follows directly from (l83i) . and by 
Fatou's Lemma 



r 

|u|| djj{u.) < liminf / ||u|| djj,°'"{u) 
Jv 



< 



The proof that Condition 2 of Definition [3] is valid follows directly from 
Lemma [Hand Lemma [51 Indeed, let us assume first that \1' G Tl C T. Then, 
by Lemma [U the functions Fj^ are weakly continuous in V, and uniformly 
bounded for u G supp/i"", a G [0,ao]- Therefore, 

lim / i^°(u)rf/i""(u) = / i^°(u)rf/i(u). (85) 
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Now, because the functions F"" satisfy the conditions of Lemma we have 

0= hm / {{I + alA)^'{u)J + uAu- B{u,u))dii'^"{u) 

= \im( [ (Ff" - F°)(u) + - FO)(u) + {Ft - F^){n)dfi--{n) 
\Jv 

+ y|^F°(u) +F°(u) + F3°(u)d/.-(u)) = J^F^iu)+F^iu) + F^in)d^^iu) 

= / {^'{u)J + uAu- B{u,u))v,v'di^{u). 
Jv 

(86) 

Now, we are ready to remove the extra hypothesis that assumed \1/ G 
Ti G T. Let \& G T be defined as in f lS^ with the function ■?/', and the vector 
fields gi G instead of ip and gj. Because V is dense in V, we have that 
for each gj, we can find a sequence |gj*''| ^ V, so that gj*^ gj, strongly 

in V, as j — * 00. Now, let us define \E'j G Tl C T by simply replacing the 
vector fields gj in the definition of ^ by gf \ i.e. 



(87) 



We also define {F^)j as in ([64]), ([65D and ([66]), with ^ replaced by ^j. 
Notice that ( |86l) holds for \l'j. It is easy to check that for every u G V fixed, 
the following convergence holds: 



which implies 



\&'(u)j ^'(u) strongly in V, as j ^ cx), 



(F;°),(u) ^ F;°(u) Vu G V^, as J ^ 00. 



(88) 



(89) 



Moreover, it is easy to see that the sequence 
there exists C > such that 



r(0 



can be chosen so that 



{F?h + (^2). + (^3°). <C||vl/'(u)ll(|f| + H|u|| + ||uf), 



uniformly in j G Now, because the right-hand side of the above ex- 
pression is integrable, by ([3]), we can use flS^ and the Lebesgue Dominated 
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Convergence Theorem to take the hmit j — oo in 



yielding 



= hm / (F^), + (F^), + (F^), d^^in) 



(^'(u), / + uAu - B{u, u))v,v' dn{u). 



V 



(90) 



This proves that Condition 2 of Definition [3] holds for every \E' G T. 
Now, we prove that Condition 3 holds. Let us define 

T = {ueV;Ei< \uf < E2} , 

and ^Z-^™) e C°°(M") as 

m 

ip^"'\ai,...,am) :=^|ajf. 

i=i 

Let us also consider a real- valued function x ^ C°°(]R), and let {w^} C 
{C°°{Q))^ n y be the orthonormal basis of H composed of the eigenvectors 
of the Stokes operator, A. Thus, for fixed m, we define for u G -B^q 



It is easy to see that \['^ G T, and that 



(u, w^))). 



(^^u) = x'E|(u,w,)r) J](u,w,)w, 



(91) 



Now, because x ^ C°°(M), we have by the Parseval identity that for every 
u, h G -ff , the following limit holds 



lim ((^^)'(u),h) = lim 



X'(j]|(u,w,)r)5^(u,w,)(h,w,) 

i=i i=i 



X'(|ur)(u,h). 
(92) 



Notice that it follows from the definition of weak convergence that supp C 
Bag, which is bounded in H. Therefore, because f G i?, we have that 
((\E'^)'(u), f) is uniformly bounded for all u G supp /i, and m G Z+. Thus, 
we may apply the Lebesgue dominated convergence theorem to obtain 

lim / ((vl/x)'(u),f)rf^(u)= / x'(|up)(u,f)d/i(u). (93) 

m^OoJy Jy 
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Now, let us denote by 1[Ei,E2] : IR — ^ M the characteristic step function on 
the interval of real numbers [Ei, E2] C M, i.e. 

1 if X e [Ei,E2l 



MeiM^) - I otherwise. 

Consider a sequence of smooth real- valued functions Xi ^ C°°(M), so that 
X'iix) l[E,,E2]ix) pointwise, as i 00, with \Xiiy)\ + \Xiiy)\ < C, for every 
i G Z+. Thus, we may again apply the Lebesgue dominated convergence 
theorem to deduce that 

hm hm / ((vl/^0'(u),f)rf/iH= /(u,f)d/i(u). (94) 

I— >oo m— >oo Jy jp 

Now, we investigate the limits for the E2 term. We obtain for every u e 
supp fi C V, 

m m 

(((vl/^0'(u),u))=x:El(u,w,)P)$^(u,w,)((u,w,)) 

i=i 3=1 

mm mm 
j=l j=l j=X j=l 

(95) 

Therefore, because supp G V, we have 

m m 

2 



lim(((v&^0'(u),u))= hm ^^(E K"' E ^^K"' = ^^d"!') II" 

(96) 

for every u G supp yU. Hence, by Fatou's lemma: 

hminf \ (((vl/^0'(u),u))d/i(u) > / xUlup) ||uf rf^(u). (97) 

Thus, we may again apply the Lebesgue dominated convergence theorem to 
deduce that 

limliminf / (((^^^'(u), u))rf/i(u) > / ||ufc//i(u). (98) 

Now, we calculate the limits for the term. Now, from ([5]), for u G supp /i C 
V^, we have -B(u, u) G V ^ and y4~^/^i?(u, u) G H. It is easy to prove that 
the following identity holds 

(5(u,u),w,)v..,y = (A-i/25(u,u), A^/V,) = (A-i/2i?(u,u),w,), (99) 
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for every j = 1,2,.... We refer the reader to [3], [U], or for a discussion 
on the properties of the negative powers of the Stokes operator. Therefore, 

(5(u,u),(vl/^0'(u))v",y 

m m 

= XiC^\{u,Wj)\'^)^{u,Wj){B{u,u),Wj)v',v 

m m 

= I (u, w,) n J](u, w,)xf{A-'/'B{u, u), w,) 



m 



m 



(100) 



I (u, w,) n j](u, av2w,)(a-v25(u, u), w,) 



= x:Ei(u,w,)n5^(Avvw,)(A-^/2i?(u,u),w,). 

i=i j=i 

Because for every u G supp /i C we have A^/^u G i/, and A^^^'^B(u, u) G 
iJ, we obtain 



hm (5(u,u),(vl/x.y(u))^,, 



hm 



X:($^|(u,w,)P)5^(AV2u,w,)(A-V25(u,u),w,) 

= (A^/^u, A-^/^B{u, u)) = (5(u, u), u)v',v = 0, 
for every u G supp /i. Therefore, by Fatou's lemma: 

hminf / (B(u, u), (^^^)'(u))y/ yc^/wlu) 

> / hminf(S(u,u),(^^0'H)y'vc^/i(u) = 0. 

m->oo 

Therefore, using ([90]), ([94]), ([98]), and ffT02D . we deduce that 
|u||^(i/i(u) < 

hm hminf / (((vI/^^'H, u)) + ((vI/^^'H, ^(u, u))v,y'rfMu) 
= hm hm / ((vl/x.)'(u),f)rf^(u)= /(f,u)d/i(u). 

r— >oo m— >oo jp 
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(101) 



(102) 



(103) 



This proves that Condition 3 of Definition [3] is vahd, concluding the proof 
that /X is indeed a strong stationary statistical solution of the Navier-Stokes 
equations. □ 



|2 . 
U IS 



lim 

n— »oo 



Now, it follows immediately from Lemma [3] that the function u 
weakly continuous in V, and bounded for u G supp for a G [0, ao]. Thus, 
we can prove strong convergence for the kinetic energy. 

Theorem 5 Given a sequence of invariant measures of the 3D Navier-Stokes- 
Voigt model, {/i""}, with a„ — *• 0, as n ^ oo, there exists a subsequence, 
denoted also by {/x""}, and a Borel probability measure fi on V, such that 

/ |u|^c//i""(u) = / |u|^c;^(u). (104) 
Jv Jv 

Now, we prove our result concerning the convergence of the net energy 

transfer. 

Theorem 6 Given a sequence of invariant measures of the 3D Navier-Stokes- 
Voigt model, {/i""}, with a„ ^ 0, as n ^ oo, there exists a subsequence, 
denoted also by {/i""}, and a Borel probability measure fi on V, such that for 
every finite pair of wavenumbers k' < k," with ^k'^k" = 0, the convergence of 
the net rate of energy transfer between k' and k," : 

lim((e«.(u))"" - (e..(u))"") = / (e.,(u) - e,.(u)) rf/x(u) (105) 

n^O Jy 

holds. 

1 1 2 

Proof. Let ao > 0. Because the function u i-^ ||uk',k"|| is weakly continuous 
in V and uniformly bounded for u G Ba^, and because a„ G [0,ao]5 for n 
sufficiently large, we can use the convergence Theorem H] to deduce 

lim ((e,,(u))°" - (e,.(u))"") = lim z/(||u,,,«.f 

n— >oo n— >oo (106) 

= u{\\u^>^^»\f) = ((eK'(u)) - (eK//(u))). □ 
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